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Pest RACT 


During gas metal are welding the metal within the heat affected zone 
undergoes a Significant thermal transient which can result in changes of 
the microstructure of the material. Such changes can significantly affect 
the strength, ductility and resistance to corrosion of the material in 
this area. The object of this study was to develop and verify a model of 
the gas metal arc process which would predict the thermal transient such 
that the mechanical and metallurgical properties of the final base metal 
a be predicted. 


The equation of heat flow in two dimensions was modified by inclusion 
of terms to account for radiation and convection heat loss. In addition, 
the fact that material properties are a function of temperature was in- 
cluded. The resulting equation was highly nonlinear, but successfully 
approximated through use of finite difference approximations. The interior 
boundary of the temperature field was the boundary of the molten weld for 
each set of welding conditions. An empirical correlation was developed, 
therefore, to predict the size of the weld pool as a function of welding 
conditions. The model was verified by actually measuring the temperatures 
with thermocouples during welding. 
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The predicted and measured temperatures were compared and found to be 
in excellent agreement whenever the weld bead was truly two dimensional. 
Peak temperatures were predicted within an accuracy of 5% or approximately 
1g2° Ff. 
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PUNT RODUETTON 


It was the objective of this thesis to develop a model of the gas 
metal are welding process such that the temperature history of the metal 
could be predicted. The gas metal arc process was chosen because it is 
a type of arc welding that is sophisticated enough to be used to join 
complex, sensitive alloys. In addition, because high deposition rates 
are possible, and automatic welding machines are available, the process 
is compatible with production welding. 

Since simple models were already in existence which work well far 
from the weld centerline, the primary area of concern was the so-called 
heat affected zone. This is the part of the base plate which, although 
it does not become molten, has been heated to a high enough temperature 
to undergo eee detectable structural change. The eee goal is 
to be able to predict the temperature history and therefore the cool-down 
rate and the final microstructure in the heat affected zone as a function 
of welding parameters. This information would be extremely useful in de- 
velopment of welding processes for new high strength steels such as 
mr-130. 

Development of the model was similar to that detailed in references 
5 and 7. First a pool size and shape correlation was developed by blowing 
out the weld puddles with high pressure argon and measuring the resulting 
holes. This information was plotted as a function of nondimensional 
parameters and an empirical relation derived. Second, the actual temper- 
atures were measured during welding operations. These values were compared 


with the results of the computer program. The input to this program 





om 
includes welding conditions and the size of the weld pool. The experi- 
mental and computed results were compared and found to be in reasonable 


agreement. 
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It. BACKGROUND 


The first real success in analytical prediction of temperature 
histories after welding was made by Rosenthal in the late 1930's. He ob- 
tained solutions for the governing heat flow equation for situations 
where the source is modeled as a point, line or plane in a quasi-steady 
state condition. Quasi-steady state here implies that the temperature 
distribution is measured relative to a moving coordinate system which 
has its origin directly under the electrode. His solutions were also 
subject to the following assumptions: 

(1) The physical and thermal properties of the medium are constant. 

(2) The radiative heat losses from the surface to the surroundings 

are small. 

(3) Any change of phase in the base metal and weld metal can be 

ignored. 

(4) As was mentioned, and most important of all, the heat source is 

modeled as a point, line or plane source. 

Experimental tests by Rosenthal and Schmerber in 1938 verified that 
the solutions found by Rosenthal were good sufficiently far from the weld 
centerline but broke down close to the weld since the molten puddle cer- 
tainly has dimensions. Rosenthal's work is summarized in reference 6. 

In the period between 1940 and 1967 many investigators tried with 
varying degrees of success to improve on Rosenthal's work. In 1952 
Wells (8) used the Rosenthal equations to predict the width of the weld 
bead. He used experimental work done by Jackson and Shrubstall to verify 


his conclusions. In 1955 Apps and Milner (2) did some work in confirmation 





-ll1- 
of the work done by Wells, but even more important for purposes of this 
work, they first expressed the theory that the total heat input to the 
base plate is composed of heat from the electric are and heat from gase- 
ous conduction from the arc column. They were also the first to use high 
pressure gas to blow out the weld pools in an attempt to prove or disprove 
theoretical size and shape predictions. Wide variance between theory and 
experimental shapes was found. Adams (1) developed engineering relation- 
ships which reasonably accurately predict centerline cooling rate and peak 
temperatures in 1958. 

One of the major obstacles to the Rosenthal solutions is that all re- 
quire the heat input to the metal to be known. It is, of course, easy to 
measure voltage and current at the weiding power supply, but, in addition, 
a value for the so-called welding efficiency must be estimated. In 1965 
Christensen et al. (3) in a detailed work measured arc efficiency. Un- 
fortunately the range of values for each welding process is very large. 

In addition in this work methods for finding length, width, depth, and 
cross sectional area of the weld bead were proposed. Christensen's sug- 
gestions for nondimensional welding parameters are also important to this 
work. 

The importance of knowing the exact temperature histories in the heat 
affected zone was summarized nicely in the 1958 Adams lecture to the 
American Welding Society by E.F. Nippes. (4) Significant changes can oc- 
cur in the microstructure in this area. As metals become increasingly 
sophisticated and sensitive to high heat these changes become even more 


important. For very high strength materials welding processes must use 
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the optimal heat input to do the least amount of damage if the strength 
of the base plate is to be preserved. For these reasons the temperature 
history in the heat affected zone must be known in detail so that the 
cool-down rate can be computed and microstructure predicted through use 
of a continuous cooling transformation diagram. This idea was suggested 
by Nippes in 1958. 

The failings of all previous models have been that they are primarily 
based on the Rosenthal solution and its inherent approximations. As was 
mentioned previously these models are good far enough from the weld center- 
line, but become inaccurate in the heat affected zone which is the area of 
primary interest. It has been proven that to make accurate predictions 
the following must be considered: 

(1) The heat source in the metal is of finite size. 

(2) The effect of the variations of thermal properties of the 

material; 

(3) radiative and convective heat loss; 

GO eediseributed teat. 

Distributed heat is that part of the energy which is transferred to the 
base plate through the gas cloud. The inclusion of these factors makes 
the original partial differential equation highly nonlinear. It has 
finally become possible to approximate the solution to such problems 
through the use of a large digital computer. 

The first truly accurate prediction of temperature in the heat af- 
fected zone was published in 1967 and 1968 by Tanbakuchi (7) and Pavelic. 


(S) Their method derives an empirical correlation for molten puddle size 





fice 
and shape, and then uses the boundary of the molten pool as one boundary 
condition in a finite difference solution. Errors of less than 10% were 
achieved. The reason for their marked success was that not only were the 
factors listed above included, but their work was limited to the gas 
tungsten are process, thus eliminating for the first time the dependence 
on an estimated welding machine efficiency. 

The principles developed by Pavelic and Tanbakuchi form the basis 
for this work. The gas metal are process was chosen for previously men- 
tioned reasons. The spray transfer mode was chosen because it 1S com- 
patible with the higher welding speeds and with more different types of 
metals. Two shielding gases were tried, argon and argon with 5% oxygen. 
The latter was chosen to improve the stability of the are and to reduce 
weld spatter. Low carbon steel was used because its properties are well 
known and it is readily available. The computer model was restricted to 
two dimensionality to reduce the computer memory required, which is very 


large even for the two-dimensional case. 
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Pit PO@eeeliZE AND SHAPE CORRELATION 


Since any sort of finite difference algorithm requires accurate 
boundary conditions this aspect of the project was very important, and 
had as its goal the derivation of empirical relationships which would 
predict points on the melting isotherm of the weld puddle. Then, since 
the melting temperature of the weld and the ambient temperature are 
known, the required number of boundary conditions are known. A theoreti- 
cal analysis is stated in some detail in reference 7. The conclusion 
reached in his work was that an analytical method of determining the 
points of maximum length and width is possible. However, this analytic 
solution is tied to the 2D solution of the heat conduction equation and 
is not confirmed by experiment. As a result, in an attempt to profit from 
previous experience detailed in reference 5, the blow-out Teed of deter- 


mining size and shape of the puddle was tried first. 


Ae Experimental Procedure 

All experimental work was done on a LINDE automatic gas metal arc 
system fed by a LINDE type SVI 300 welding power supply. After considera- 
ble trial and error a satisfactory blow-out system was constructed and is 
depicted schematically in Figure 1. The copper tubing was attached to the 
torch weer welding machine so that it had a constant orientation with 
respect to the are and the molten puddle. The flow restricting valve was 
a coarse regulating needle valve and was left partially open. The blowing 
out of the puddle was accomplished by stopping the welding machine and 
opening the gate valve simultaneously. This apparatus was simple and uwn- 


sophisticated but was found to be effective once sufficient experience 
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running the welding machine had been accumulated and once the correct 
angle of the final length of copper tubing was determined. The pressure 
of the argon used to blow the puddles out was not measured but is esti- 
mated to be 500 pounds per square inch. 

All the blow-out specimens measured two inches by 10.8 inches. This 
size specimen was chosen because it was narrow enough to be almost free 
of distortion during the welding anc yet wide enough to prevent heat re- 
flected at the edges from affecting the welding conditions. Experience 
also showed the particular welding machine used consistently produced a 
Stable weld after an initial run of three to four inches. Therefore a 
length of over ten inches allowed some choice of the exact moment when 
the puddle was blown out. 

The following table summarizes the welding .conditions for the blow- 


out specimens. 


Voltage 15-20 volts 

Current 190-280 amps 

Speed 28.5 to 48.7 inches per minute 
Thickness Wy USselg ako vele 

Material mild steel 


The driving parameter is the speed. It was found that at each new ‘speed 
setting the voltage and current must be adjusted very exactly to get a 
two-dimensional puddle and yet to prevent the puddle from dripping through. 
Attempts were also made to get two-dimensional puddles in 1/8th inch steel, 


but the increased width of the weld puddle allowed the molten material to 





Sve 
drip through before it could be blown out. This result is in agreement 


with Pavelic's. Experimental results are included in Appendix C. 


Be Empirical Correlation 
By manipulation of the line source solution to the heat conduction 
equation a nondimensional welding parameter can be derived and is de- 


mined as follows: 


ee oe EL 
27QtH,, 


Details are included in Appendix D. This parameter is a function of the 
Material properties and of the welding parameters of primary interest, 
namely voltage and current. It is not so obvious that the above parameter 
provides a method of excluding welding efficiency from further considera- 
tion for this procedure. If the dimensions of the puddle can now be ex- 
pressed in terms of N*, then for a particular machine, the isotherm of 
melting can be determined and the remainder of the solution is a standard 
finite element procedure consisting of averaging the temperatures of sur- 
rounding nodes until the solution converges. 

There is one flaw in this reasoning, however, and that is that the 
effect of the welding speed has not been included. A coordinate system 


is defined as follows: 
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where w is the moving coordinate in the direction of the welding and y 
is the coordinate used to express distance from the weld centerline. 
These coordinates can be nondimensionalized through dimensional analysis 


as follows: 


* Vy 


we = )Vw 


< 
i 


and the effect of welding speed has been included. 

In a manner similar to that used by Pavelic the blown-out weld pud- 
dles were measured using a Nikon machinist's microscope. These results 
were fed into the service program "puddle fit" which is included in Ap- 
pendix B. As can be seen from Figures 2 and 3 the data falls roughly 


into an area which might be described by a relationship of the form 





30 


10 


e Good 


x Fair 


o Poor 


Figure 2. 
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Figure 3. MAXIMUM PUDDLE WIDTH CORRELATION 
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where C is the intercept of a straight line on a log log plot and m is 
its slope. Program "puddle fit" performed this operation using a least 
squares approximation. The results of this program are the following 
relations: 


for maximum length 
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ee eo ee COMPUTER MODEL 


This section is a review of the computer model and a discussion of 
the minor changes made to it. The basic equations are unchanged. That 
is the model retains correction terms for diffused heat, convection, 
radiation, and variable material properties. Changes were of a minor 
sort and were intended to, first of all, make the program run in a speed 
range higher than that considered previously and second, to reduce the 


actual run time of the entire program. 


A. Basic Equations Used 

If it is assumed there is no heat generation in an isotropic solid 
and that the thermal properties are a function of temperature, the dif- 
ferential equation of heat may be written 

Ci) te (4-1) 
Ds 

for a standard Cartesian coordinate system. If a quasi-stationary state 
is defined in which temperature is not a function of time and the coordi- 
nate system is measured relative to the moving are this equation is 


changed to the following: 


Ve: (kT) 2-peV— (4-2) 
W 
where w is the coordinate in the direction of motion and is described by 
the following: 
w= x - Vs 
If a normalizing factor is defined 


x, = = “min | Gea) 
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where S is a constant and will be discussed in detail later. Then all 
dimensions can be normalized by dividing by Xo Similarly temperatures 
are normalized by dividing by the difference between melting temperature 


and ambient temperature, that is, 


: to LO 
Ts = 
Tm -To 


After normalization, rearrangement, and introduction of the two- 


dimensional approximation equation (4-2) becomes: 


?, ¢, 9, ?, 2 ?, 
g2 re 92T Vx aT 1 ak (T*) Ba Bal 








' 
© 


(4-4) 


4 t pa CSRS ’ Tae a 
Qwe2 8 dy? TH) Ow =k (TH) OTH awit ay 


However, aS waS clearly shown by the differences in reference 5 and 
7, not all the energy is transferred into the base plate through the molten 
pool. Some energy is transferred directly from the ionized shielding gas 
to the base plate. A function which approximates this distributed heat is 
given below: 
an(r) = q (0)e°™ 


2 a 


where r2 =wt ty 


and which when normalized becomes 


q(0) 


-Cx,2 (wit? + y#2) 
(Tm = Wey) =) 


dn ws, y*) = 
Losses due to convection and radiation can also be included and the result 
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In accordance with standard procedures the partial differentials are 


now replaced by finite differences which are now defined. 


oT* OT eaeea - a@i-1,j) _x ai-1,5 ‘i: 
rls ee : Ny ay re .- 
a Sie jci 2d ai,j (44,jt4i-1,5) “itl.3 
Gita re 


ee —=" 
eee, 42 22s 





ee ei of oe 
eee PRG PPP 5) 3.942 
ies nt 
}- oer aoe Tye 
b; 5) Jk (bi 5 ti bs 5-1? >) ilk 
ger 9 ’ 5 ‘ 
awez * oS od, ae (Ce area a ee 
“od iti OR te ae eg 
2, wv 
es ° 
Cay ih 


f, 


~ oo 2 * 
in Sanat eGo ee) us x 
= ne Pe ail bi (Di 5 id Pa ial ~2 Jt 


2 v's 
pa. . = ee 
ae = 15) 1,J-l > J 


where the computational module is defined in Figure 4. 
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Figure 4. 


Substitution of these expressions into equations and rearranging 


yields the following: 
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B. Convergence and Uniqueness 
In reference 7 Tanbakuchi used the following criteria for grid spac- 


ing parameters in order to insure convergence: 





i 20min 
<--_- 
max om 
ie h Which is a nondimensional grid spacing parameter in the program 


Max 


itself, is less than twice the minimum value of diffusivity divided by 
another value of diffusivity, then the solution will converge and will 


be unique. If Ree is arbitrarily set to 1, then 


a = 2a, 
‘o) min 


and convergence and uniqueness will be satisfied if the length normalizing 


meetor 1s defined 


- 2onin 
Se (4-7) 





However, it was found that in the speed range greater than approxi- 
mately 30 inches per minute the lenzth normalizing factor became very 
small which caused the puddle dimension to become very large. As a result, 
computer limits for division by very small numbers were exceeded in sub- 
routines TYPFN or COLUMN because Ko bessel function values returned from 


subroutine BESSEL were very small. Similar problems had to be corrected 





Be a 

throughout the program as a result of the increased welding speed. The 
solution to this problem was to increase the length normalizing factor, 
X52 by replacing the 2 in equation (4-7) by an arbitrary factor (see 
equation 4-3). This factor can be set at a value higher than 2 if the 
grid spacing parameter is set proportionally lower than 2. For most 


runs S equaled 4 and h equaled .5. 


max 
Unfortunately, however, such careful work does not insure convergence. 

As is pointed out by Pavelic, equation (4-4) is highly nonlinear and there- 

fore a finite difference approximation need never converge. The experi- 


mental results are available, however, for ready comparison so that a 


diverging solution should be immediately obvious. 


on Changes in the Computer Program 

As was mentioned at the beginning of this section, the changes made 
to the program were minor, but worth mentioning. The first change was in 
the definition of the normalization factor and has been discussed in some 
detail. The second change was the establishment of two holding arrays in 
Which successive iterations of the temperature surface are stored. After 
each complete iteration each element in each array is compared with the 
corresponding element in the other. If every element is within a toler- 
ance which can be specified in the input deck, execution stops and results 
are printed out. This improvement in the program reduced run time on an 
IBM 360 model 155 from several minutes to approximately one-and-one-half 
minutes. However, all the changes made so little difference to the struc- 


ture of the program that the flow charts in reference 7 still apply. 
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V. EXPERIMENTAL MEASUREMENT OF TEMPERATURE HISTORY 


Experimental temperatures were measured using the same welding ap- 
paratus described in section 3. Specimens were either 1/16th inch or 
1/8th inch cold rolled SAE 1018 steel and measured 17.5 inches long by 
eight inches wide. Six thermocouples made of 24-gauge chromel and alumel 
wire were placed in each specimen. Filler metal was Airco type 675 which 
is recommended for low-carbon steel, and argon plus five percent oxygen 
continued to be used as the shielding gas. A Honeywell type 1508 aeeoree 
er was used to measure and present the output of the thermocouples. A 


schematic of the measurement circult is presented in Figure 5. 


A. Thermocouple Emplacement 

Welding conditions were first established using scrap steel. Then 
program puddle fit was used to estimate the width of the weld bead. Then 
thermocouples were placed at varying distances from the weld centerline, 
starting just outside the predicted width of the weld bead. Since temper- 
atures in the base plate can be predicted accurately by the point source 
theory as close as 1/4 of an inch from the weld centerline, all thermo- 
couples on all three temperature measuring runs were placed at distances 
less than .25 inches from the weld centerline. 

Since 24-gauge wire 1S approximately .020 inches in diameter, thermo- 
couple holes were drilled with a diameter of approximately .040 inches. 
After the holes were drilled the wire was faced through the plate and the 
ends of the two wires were melted together. The small bead thus formed 
was hammered into the top surface of the plate to insure an intimate con- 


tact with the entire surface of the hole. After the thermocouples were in 
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Dlace the actual distance from the weld centerline to the centerline was 
measured using the microscope mentioned previously. Figure 6 depicts a 


typical thermocouple installation. 


B. Temperature Measurement Procedure 


Figure 7 depicts the welding machine during a temperature measurement 
run. Voltage and current were measured by meters on the power supply. 
Speed was measured by using a stop watch to time travel over a known dis- 
tance. To start a run the Honeywell recorder was started at a speed of 
one inch per second. It is a feature of this recorder that timing marks 
are automatically made on the paper. Next the welding machine was started, 
and the start point of the welding machine was noted on the recorder paper. 
The welding machine was stopped several inches after the last thermocouple 
was passed, and the recorder after all temperatures had decayed sufficient- 
ly. In addition, for program "experimental data," the distance from the 
start of the weld to the first thermocouple and the distances between 
thermocouples were measured. 

Three experimental runs were made with varying degrees of success. 


Welding parameters for each are listed below in Table l. 


TABLE 1. 
Run al Z 3 
Voltage (volts) S EGS 22 
Amperage (amps) 190 235 CW 
Thickness (in.) 1/16 7 M6 ye 


Speed (ipm) 28.5 43.5 28.8 
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See 
oe Analysis of Experimental Data 
The output of the Honeywell recorder is in the form of deflections 
continuously displayed. A sample curve is reproduced in Figure 8. The 
deflection was measured every half second and converted to a temperature 
by use of a calibration curve. These temperatures for an experimental 
run are included in Appendix C. This data plus the distances already 
mentioned form the input to program "experimental data" which puts the 
experimental data in the same coordinate system as the computer output 
for comparison. A listing of program experimental data is contained in 


Appendix B. 
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Poy, COMPA RSON OF COMPUTER AND EXPERIMENTAL RESULTS 


Welding parameters for the three experimental runs are listed in 
Table 1 in section V. Figures nine through 15 are a comparison of the 
temperatures measured in runs 1 and 3 with those computed using the weld- 
ing parameters and the thermocouple location. For information the thermo- 
couple locations are included in Appendix C. 

If error is defined: 


io dhe 


R Se 100% 
ERROR 6 ( Z ) 


then the error ranges in magnitude from as little as five percent at the 
maximum separation for thermocouple four to as great as 18 percent at 

the maximum separation for thermocouple one. This represents a differ- 
ence in temperature ranging between 132° R. and 476° R. at the points of 
maximum separation. With this particular definition of error, the error 
observed is somewhat higher than that claimed by Pavelic in his work. It 
must be remembered, however, that his definition of error is based on a 
comparison of the peak temperatures regardless of when they occur rela- 
tive to each other. Using his definition the magnitude of the observed 
error is decreased considerably. 

As Pavelic pointed out the ideal comparison between experimental and 
computed data would occur if the curves were to overlap. However, in 
practice this is difficult to accomplish. Since experimental results 
generally lag the computed results, the contact between the wall of the 


hole and the thermocouple is immediately suspect. In fact, however, the 
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problem may be in the numerous factors and variables which are a part of 
welding and yet are extremely hard to quantify. 

As was mentioned before, welding conditions were determined first by 
choice of speed. The speed forced an optimal choice of voltage and cur- 
rent since the desired goal was a two-dimensional weld bead. However, 
as the speed or thickness increased, it was observed that the surface 
of the molten puddle flattened out and the resulting bead departed rad- 
ically from true two dimensionality. Experimental runs one and two were 
sectioned and macroscopically etched to determine the shape of the weld 
bead cross section. The results are represented in Figure 16. Since the 
thermocouples were a discontinuation in the specimen surface, all thermo- 
couples came in contact with the molten metal and consequently all gave 
relatively high readings. Figure 15 is thermocouple three of run three.. 
Using either definition, the error for this thermocouple is approximately 
20 percent and is unacceptably high, and this was the best thermocouple 
mem this run. 

The problem in these two runs is that too much power is being trans- 
ferred to the metal. All welding handbooks consulted advised reducing the 
current and the voltage to reduce spatter and to reduce the tendency of 
the weld bead to spread out at the top. Of course, if the power were re- 
duced, the penetration is reduced and two dimensionality is lost anyway. 
One other possible variable is the type of welding gas. Argon plus 5% 
oxygen was used, however, which is recommended for gas metal arc welding 
on low-carbon steel. 


The conelusion reached was that the method first used by Tanbakuchi 
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and Pavelic is applicable to gas metal arc welding in a spray transfer 
mode. However, a lesser degree of success was achieved at very high weld- 


ing speeds. 
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TWO-DIMENSIONAL WELD BEAD (SLOW SPEED) 


NON-TWO-DIMENSIONAL WELD BEAD (HIGH SPEED) 


Figure 16. 
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VII. SUMMARY AND RECOMMENDATIONS 


This method of prediction of temperature history in the heat af- 
fected zone is certainly useful because it 1S accurate. However, be- 
cause of the constraint of two dimensionality it is, in the opinion of 
the author, of only limited usefulness. For any sort of temperature pre- 
diction process to have any practical usefulness it must predict the 
temperatures under the actual welding conditions. Actual welding condi- 
tions in this case imply a somewhat lesser heat input per unit time than 
is necessary to achieve two dimensionality. The obvious implication is 
that this method should be extended into three dimensions. Only specula- 
tion is possible, however, whether the practical aspects of this method 
could be applied in three dimensions. It is doubtful, for example, that 
the blow-out system would work as well as it did. Since the blow-out 
equipment is obligated to be angled down toward the weld pool, it might 
be impossible to determine the depth of the weld pool to the high degree 
of accuracy required. 

The obvious answer to these difficulties is to make the angle between 
the plate and the blow-out nozzle as small as possible and increase the 
pressure of the blow-out gas, but only actual attempts would determine if 
the method is feasible. 

Unfortunately also, even the two-dimensional version used for this 
analysis required 400,000 bytes of memory to analyze an area one inch by 
three inches. It is obvious that this requirement grows rapidly as either 
more area is considered or as the third dimension is included. The need 


to have a large sophisticated computer greatly restricts the practical 





Sige 
value of this analysis since most production facilities like shipyards 

or aircraft manufacturing plants do not have such an installation. Of 
course if peripheral devices are available, they could be used to augment 
in core memory, but only if the program went through a major revision. 

With the advantage of hindsight, a far more reasonable approach seems 
to be to model the molten pool as a shape composed of many point sources. 
These point sources would be analytically varied in magnitude through the 
value of the molten metal to give the proper shape and heat input to the 
metal. This idea is similar to the analytical sources and sinks used to 
model shapes in the potential flow of fluids. The advantage here, if 
such an idea is possible, is that the temperature could be computed at any 
point in the material without an enormous amount of computer memory being 
required since the temperature at every point is simply a combination of 
the temperatures resulting from the action of each small source. It is 
obvious that the above statement assumes linearity which would have to be 
checked experimentally. However, for practical applications, it is the 
opinion of the author that a way other than finite difference solution must 


be found in order to reduce the computer memory requirements. 
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B< Experimental Temperatures 


Run #1 
Distance from 
Local Origin (in.) Temperature (°F.) at Thermocouple 
if 2 3 Ly 5 6 
O50 17 Ae Ta ad qT 77 
Ce) 50 330 E20 90 92 90 
lee @) 875 730 240 205 290 170 
DS JEL Se OCS 610 490 605 ooo 
ZeQ WLS TANGO 810 690 25 605 
ae 1370 i769 960 790 840 MS, 
3.0 1410 Lei 1065 890 920 805 
Seo 1410 376 1095 970 970 890 
4.0 1380 Abs ishs 75 JLO.IUG: 1005 955 
4.5 fess ees 1190 eee Toei, 990 
50 1330 1305 1190 O70 1045 1018 
Sar ALPS 1290 aces LOS LOSe 1030 
Grail) 1270 1269 ili bei: 1085 1050 1040 
er) 1240 1240 LG 5 170 1040 1050 
7.0 126 i205 1145 1065 1030 LOSS 
es: 1190 1180 LSS ESO O25 1050 
B20 algae: JAC Zs 1040 1020 1040 
Beno nh Mesye: 1130 OS LO ZS 1010 1040 
9.0 a NEN: RES 1035 1010 995 1030 
9.5 1090 1090 1080 1000 990 1025 
Distance from 
Weld C/L. (in.) plead) ~139 <S5 Suilyoial . 190 2220 
Meco RTC i) 9.74 2 re 78 ISAs 2s 24.95 


XCCHRT(1) is the distance in inches from the origin of the run to 
the local origin of each curve. 
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The Rosenthal solution for the 2-D case assuming a thin plate and a 
line source is: 


il AVw - 


as Ka) (D-1) 


on ee 
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where Q = 3.415 nEI (Btu/hr) 
and r = vw2 + y2 (ft.) 


Hn? the heat content of the material, is defined: 


Hee cp (T= T 


m 0) 


which ignores the heat required for the pearlite-austenite transformation 
and the latent heat of fusion which are small effects and very difficult 
to measure accurately. By multiplying and dividing equation (D-1) by 
Oc oe - wey the following results: 

te =a = 


2. SLE 1 lee dog m~ ~0 : 7 
4 a Se oe eee Ko (AVP) (D-2) 
k pe(T - To) 


When reduces to: 
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when the nondimensional quantities defined in sections III and IV are in- 


cluded. If a nondimensional welding parameter is defined: 
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then the Rosenthal solution in the nondimensional form is the following: 


A 
-We 


T= oN A Car) 


The same parameter was used for this study since it combines plate 
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thermal properties, thickness, welding voltage, and welding current in a 


convenient form. 








Thesis 4a 
L675 Lipfert 134938 
Temperature histories 


in thin steel plate dur- 
ing gas metal arc we ld- 
ing. 





